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We investigate continuous time random walks with truncated Q-stable trapping times. We prove 
distributional ergodicity for a class of observables; namely, the time-averaged observables follow 
the probability density function called the Mittag-Lefiler distribution. This distributional ergodic 
behavior persists for a long time, and thus the convergence to the ordinary ergodicity is considerably 
slower than in the case in which the trapping-time distribution is given by common distributions. 
We also find a crossover from the distributional ergodic behavior to the ordinary ergodic behavior. 
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7 The ergodic theorem ensures that time averages of ob- 

8 servables converge to their ensemble averages as the av- 

9 eraging time tends to infinity. On the other hand, a dis- 

10 tributional ergodic theorem states that the probability 

11 density functions (PDFs) of time averages converge to the 

12 Mittag-Leffler (ML) distribution. This property is called 

13 infinite ergodicity in dynamical system theory, because 

14 it is associated with infinite invariant measures [l| . Fur- 

15 thermore, in recent years, the distributional ergodic prop- 

16 erty has been found for some observables in stochastic 

17 models such as continuous time random walks (CTRWs) 

18 For example, the time-averaged mean square dis- 

19 placement (TAMSD) [Eq. ^] for CTRWs is a random 

20 variable even in the long measurement time limit and its 

21 PDF follows the ML distribution. It has been pointed out 

22 that this distributional ergodic behavior is reminiscent of 

23 the observations in biological experiments that showed 

24 that TAMSDs of macromolecules are widely distributed 

25 depending on trajectories 0]. In addition to these bio- 

26 logical systems, CTRW-type systems are used to explain 

27 a broad range of phenomena such as charge carrier trans- 

28 port in amorphous materials Q Jracer particle diffusion 

29 in an array of convection rolls j5[, and human mobility 

30 0. 

31 One of the important problems on stochastic models 

32 such as CTRWs is to clarify the condition of the dis- 

33 tributional ergodicity. It has already been known that a 

34 few observables including the TAMSD show the distribu- 

35 tional ergodicity in CTRWs. But any general criterion 

36 for an observable to satisfy the distributional ergodicity is 

37 still unknown. Another important problem to elucidate is 

38 finite size effects For CTRW-type systems, a power 

39 law trapping-time distribution is usually assumed, and 

40 thus rare events — long-time trappings — characterize the 

41 long-time behavior. These rare events, however, arc often 

42 limited by finite size effects. For example, if the random 

43 trappings are caused by an energetic effect in complex 

44 energy landscapes, the most stable state has the longest 

45 trapping time, thereby causing a cutoff in the trapping- 
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46 time distribution. In fact, for the case of macromolecules 

47 in cells, the origin of trappings is considered to be ener- 

48 getic disorder: strong bindings to the target site, weak 

49 bindings to non-specific sites, and intermediate bindings 

50 to sites that are similar to the target site Q. Because 

51 the binding to the target site should be most stable with 

52 the longest trapping time, there must be a cutoff 

53 Similarly, if the trappings are due to an entropic effect 

54 such as diffusion in inner degrees of freedom (diffusion 

55 on comb- like structures is a simple example @; see also 

56 [l^l), the finiteness of the phase space of inner degrees 

57 of freedom results in a cutoff. The CTRWs with such 

58 a trapping-time cutoff show distributional ergodic fea- 

59 tures for short-time measurements, and become ergodic 

60 in the ordinary sense for long-time measurements. But 

61 this transition from distributional ergodic regime to or- 

62 dinary ergodic regime has not been elucidated. 

63 In this study we employ a truncated one-sided stable 

64 distribution as the trapping-time distribution, and 

65 show that the distributional ergodic behavior persists for 

66 a remarkably long time compared to the case of common 

67 distributions with the same mean trapping time. We also 

68 show that the time-averaged quantities for a large class of 

69 observables exhibit the distributional ergodicity. As an 

70 example, numerical simulations for a diffusion coefficient 

71 are presented. We use the exponentially truncated stable 

72 distribution (ETSD) prop osed in and the numerical 

73 method presented in This ETSD is useful for rig- 

74 orous analysis of transient behavior, because it is an in- 

75 finitely divisible distribution [T^j and thus its convoluted 

76 distribution or characteristic function can be explicitly 

77 derived [Eqs. ^ and dH)]. 

78 Truncated one-sided stable distribution. — In this study, 

79 we investigate CTRWs on d-dimensional hypercubic lat- 

80 tices. The lattice constant is set to unity, and for simplic- 

81 ity, the jumps are allowed only to the nearest-neighbor 

82 sites without preferences. Let r{t') G Z"^ be the position 

83 of the particle at time t' . Moreover, we assume that the 

84 successive trapping times [k = 1,2, ...) between jumps 

85 are mutually independent and the trapping-time distri- 

86 bution is the ETSD Ptl (t. A) defined by the canonical 
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form of the infinitely divisible distribution [T^ ]: 

-OO 

V'(C,A) = / (e<^-l)/(T,A)dr, 



(t<0) 
(r > 0), 



The function /(t, A) is defined by [12 
0, 



/(t,A) = 
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where T{x) is the gamma function, c > is a scale factor, 
and a G (0, 1) is a constant. The parameter A > char- 
acterizes the exponential cutoff [Eq. ([6])]. When A = 0, 
Ptl(''', A) is the one-sided a-stable distribution with a 
power law tail Ptl(t,0) - 1/t1+" as t -!> cx) [li]. The 
function V'lCi ^) can be expressed as follows: 

^(C,A) = -c[(A-zC)"-A"]. (4) 

Hence, we obtain n^{(, A) = ip{n^/°'(, n^^°'X), where n > 
is an integer. Therefore, if (fc = 1, 2, ...) are mutually 
independent random variables each following Ptl (Tfe , A) , 
then the n-times convoluted PDF Pxl(t, A), i.e., the PDF 
of the summation T„ = X]fc=i '''k, is given by 
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FIG. 1. (Color online) TAMSD {Sx)2(A,t) vs time interval A 
(log-log plot) for the one-dimensional system (d = 1). Total 
measurement time t is set as t — 10'' and other parameters as 
A = 10"^ a = 0.75, and c = 1. The TAMSDs are calculated 
for 17 different realizations of trajectories; a different symbol 
corresponds to a different realization. The solid line is their 
ensemble average. 



117 for long t, where Nt is the number of jumps until time t. 

118 From this equation, we find that ht behaves similarly to 

119 Nt ■ It is important that many time-averaged observables 

120 for CTRWs can be defined by Eqs. Q and For 

121 example, the TAMSD, 



This is an important outcome of the infinite divisibility 
and makes it possible to analyze transient behavior of 
CTRWs. Moreover, from Eq. (|4]) and the inverse trans- 
form of Eq. ([T|), we obtain an explicit form of PTh{T, A) 
through the similar calculation shown in : 
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\rit' + A)-r{t')\^dt', (10) 
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Ptl(t,A) = — 



T{ka 

k=l 



kl 



— ^ (— CT sin(7rfca). 

(6) 

CTRWs with truncated a-stable trapping times. — Now, 
we consider the time average of an observable h{t'): ht = 

dt'h{t')/t, where t is the total measurement time. We 
assume that h{t') can be expressed as 



122 can be approximately obtained by the time average of 

123 h{t') with Hk defined as Hk = A+2 J2'^~l {drk-dri)9{A~ 

124 (Tk—Ti)), where d-dimensional vector drk is the displace- 

125 ment at the time Tfc, and 6{t) is defined by 9{t) = t for 

126 t > 0, otherwise 9{t) = 0. It is easy to see that (Hk) ~ A 

127 and (HkHk+n) - (Hk) {Hk+n) = for n > 1. Using 

128 Eq. we have 



(5r)2(A,t)ci AiVt/t. 



(11) 
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129 From Eq. (fTTj) . we obtain a relation between Nt and the 

130 diffusion coefficient of TAMSD as 



k=l 



where Tk > (fc = 1,2,...) is the time when the fc-th 
jump occurs, and Hk {k = 1, 2, ...) are random variables 
satisfying {Hk) = (H) and the ergodicity with respect to 
the operational time fc. 
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To satisfy Eq. (|5]), the correlation function {HkHk+n) — 
(Hk) (Hk+n) should decay more rapidly than n~'^ with 
some constant 7 > [13]. It follows from Eqs. ([7|) and 
dl) that 
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In Figll] TAMSDs calculated from 17 different trajecto- 
ries are displayed as functions of time interval A. This 
figure shows that the TAMSD grows linearly with A, 
and the diffusion coefficient Dt is distributed depending 
on the trajectories. 

PDF of time- averaged observables. — In this section, we 
derive the PDF of time-averaged observables ht. Because 
ht and Nt have the same PDF [Eq. p2|) ]. we can study 
Nt instead of ht- We have the following relations: 

G(n; t) = Prob {Nt < n) = Prob (T„ > t) 

= Prob 1^^^ Tk>t^, (13) 
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where Prob(-) is the probabihty and is the trapping 
time between (k — l)-th and fc-th jumps (fc = 1,2,...). 
Prom Eq. (jlSp . we obtain 
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175 where we have defined 
dte~^'^G{n\ t) and used 



G{n; s) as G(n; s) = 
dte-^'Gin-.t) and used Eq. (g]). Next, we de- 
fine a function g{n] s) by ^(n; s) := G{n + 1; s) — G(n; s). 
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where we have used Eq. ([5]) and the fact that {k ~ 
1,2,...) are mutually independent. Furthermore, we 
change the variables from n to z as n = with t being 
set. Then, by using Eqs. dH), (fT3|) and (fT4| . we have 




where ak{t) [k = 1,2,...) is defined by ak{t) = 

/(,^ dre'*'^'^ ^ . Differentiating Eq. ([T5|) with respect 
to 2, we have the PDF of z = Nt/t°': 

T{ka + 1] 



Diffusion coefficient D, 



FIG. 2. (Color online) The PDF of the diffusion coefficient 
Dt for d = 1. Each PDF is normalized so that its me an value 
equals unity. Dt is calculated from TAMSD (Sxy^{A,t) by 
least square fitting over the interval < A < 10. The results 
for three different values of measurement times are presented: 
t = 10* (circles), 10* (squares), and 3 x 10* (triangles). The 
other parameters are set as A = 10"'', a — 0.75, and c — 1. 
The lines correspond to the theoretical predictions given by 
Eq. (|16p . Note that no adjustable parameters were used to 
obtain these figures. 



Then, by taking a (discrete) Laplace transform with 
respect to n, X^J^o 6~"''5("-; s) = g{v, s), we have 



Because of Eq. (l9|), Eq. dm) is the PDF of the time- 
averaged observables htt^^" / (H) including the diffusion 
constant Dtt^~°' [Eq. (|T2|) ] as a special case. When A — 0, 
the PDF /o(z), which is the ML distribution Q, is time- 
independent. Namely, the time-averaged observables are 
random variables even in the limit t oo; this property 
is called the distributional ergodicity. On the other hand, 
when A > 0, the PDF tends to a delta function. Thus, the 
time-averages converge to constant values as is expected 
from the ordinary ergodicity. The PDF of Dt is shown 
in Fig. [2] for three different measurement times t. It is 
clear that the PDF becomes narrower for a longer t. The 
analytical result given by Eq. ([TS]) is also illustrated by 
the lines. 

Relative standard deviation. — Next, in order to quan- 
tify a deviation from the ordinary ergodicity, we study a 
relative standard deviation (RSD) of time-averaged ob- 
servables R{t) = {{htY)cl (ht), where (■) is the ensem- 
ble average over trajectories and (•)^ is the corresponding 
cumulant. If i?(t) w 0, the system can be considered to 
be ergodic in the ordinary sense, whereas if R{t) > 0, the 
system is not ergodic. To derive an analytical expression 
for R{t), we take the Laplace transform of Eq (IT4l) with 
respect to t: 
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where we used an approximation by assuming s,A <C 1. 
From Eq. (|18p , we can derive arbitrary order of moments 
of Nt- For example, the first moment (Nt) is given by 
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cF(a + 1) 
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t » 1/A. 



(19) 



cA"-ia 2cA"q;' 
The ensemble-averaged mean square displacement 
(EAMSD) for CTRWs is known to be proportional to 
{Nt) 0: ((5r)2)(t) ^ (Nt). Thus, the EAMSD of the 
present model shows transient subdiffusion, i.e., subdif- 
fusion for short time scales and normal diffusion for long 
timescales Similarly, the second moment can be de- 
rived and we have the RSD for Nt as follows: 



(Nt) 




1, 



t < 1/A 



2r3(a+l) _ 

r(2a+l) " . ^ V- (20) 

l-a)i/2A-i/2i-i/2^ t»l/A. 



Note that the RSDs for ht and Dt also follow the same 
relations, because they differ only in the scale factor 
[Eqs. ^ and ([T^ ]. From these results, the crossover 
time tc between the distributional and ordinary ergodic 
regimes is given by 

(1-a) 



2r2(a + l)/F(2a + l)-l^ ^^^^ 

As shown in Fig. [31 the RSD remains almost constant 
before the crossover time tc, and starts to decay rapidly 
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FIG. 3. (Color online) RSD ^JTpf)J (^*) vs total measure- 
ment time t for d = 1. Dt is calculated from the TAMSD 
((5a;)2(A, t) by least square fitting over the interval < A < 1. 
Three different values of A are used: A — 10"'* (circles), 
10~^ (squares), and 10~^ (triangles), and a and c are set as 
a — 0.75 and c — 1, respectively. The lines correspond to the 
theoretical prediction given by Eq. (|20p : the solid line is the 
result for short time scales t <^ tc, whereas the dashed lines 
are for long time scales tc- The intersections of the solid 
and dashed lines correspond to the crossover times tc given 
by Eq. (|21[l . The pluses are the RSD for the case of the ex- 
ponential distribution P{t) = exp(r/ (r))/ (r) with the same 
mean trapping time as the ETSD with A = 10^'' (triangles): 
(t) = cA"~^a. The dot-dashed line is a theoretical prediction 
for the exponential distribution: R{t) = (caA'^^^/t)^''^. 



after the crossover. In Fig. [3l the RSD for the exponen- 
tial trapping-time distribution which has the same mean 
trapping time (t) as the ETSD with A = 10~^ is also 
shown by pluses. It is clear that the RSD for exponential 



distribution (pluses) decays much more rapidly than that 
for the ETSD (triangles). 



203 Summary. — In this study, we have investigated the 

204 CTRWs with truncated a-stable trapping times. The 

205 three main results are as follows: (i) We proved the distri- 

206 butional crgodicity for short measurement times; namely, 

207 the time averages of obscrvablcs behave as random vari- 

208 ables following the ML distribution. Moreover, we de- 

209 rived the PDF at arbitrary measurement times. It is very 

210 interesting to compare this analytical formula [Eq. (1161) 1 

211 with the results for lipid granules reported recently [l6|. 

212 Wc should also note that the limit distributions, the ML 

213 distribution for the case of observables studied in this 

214 paper, depends on the definition of observables [13] • (ii) 

215 We found that the distributional ergodic behavior per- 

216 sists for a long time. In other words, the convergence 

217 to the ordinary ergodicity is remarkably slow in contrast 

218 to the case in which the trapping-time distribution is 

219 given by common distributions such as the exponential 

220 distribution. This indicates that, in real experiments, 

221 the time-averaged quantities could behave as random 

222 variables even for considerably long measurement times. 

223 (iii) We found a crossover from the distributional ergod- 

224 icity in the short-time regime to the ordinary ergoodicity 

225 in the long-time regime. Finally, it is worth mentioning 

226 that these three main results are valid for a large class 

227 of observables. This implies that it is possible to choose 

228 an observable which is easy to measure experimentally. 

229 Then, the system parameters a and A can be expcrimen- 

230 tally determined by the short- and long-time behavior 

231 of the RSD R{t) [Eq. respectively. 
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